Genomic data resources of the Brain Somatic Mosaicism Network for neuropsychiatric diseases.

Somatic mosaicism is defined as an occurrence of two or more populations of cells having genomic sequences differing at given loci in an individual who is derived from a single zygote. It is a characteristic of multicellular organisms that plays a crucial role in normal development and disease. To study the nature and extent of somatic mosaicism in autism spectrum disorder, bipolar disorder, focal cortical dysplasia, schizophrenia, and Tourette syndrome, a multi-institutional consortium called the Brain Somatic Mosaicism Network (BSMN) was formed through the National Institute of Mental Health (NIMH). In addition to genomic data of affected and neurotypical brains, the BSMN also developed and validated a best practices somatic single nucleotide variant calling workflow through the analysis of reference brain tissue. These resources, which include >400 terabytes of data from 1087 subjects, are now available to the research community via the NIMH Data Archive (NDA) and are described here.

Candida albicans selection for human commensalism results in substantial within-host diversity without decreasing fitness for invasive disease.

Candida albicans is a frequent colonizer of human mucosal surfaces as well as an opportunistic pathogen. C. albicans is remarkably versatile in its ability to colonize diverse host sites with differences in oxygen and nutrient availability, pH, immune responses, and resident microbes, among other cues. It is unclear how the genetic background of a commensal colonizing population can influence the shift to pathogenicity. Therefore, we examined 910 commensal isolates from 35 healthy donors to identify host niche-specific adaptations. We demonstrate that healthy people are reservoirs for genotypically and phenotypically diverse C. albicans strains. Using limited diversity exploitation, we identified a single nucleotide change in the uncharacterized ZMS1 transcription factor that was sufficient to drive hyper invasion into agar. We found that SC5314 was significantly different from the majority of both commensal and bloodstream isolates in its ability to induce host cell death. However, our commensal strains retained the capacity to cause disease in the Galleria model of systemic infection, including outcompeting the SC5314 reference strain during systemic competition assays. This study provides a global view of commensal strain variation and within-host strain diversity of C. albicans and suggests that selection for commensalism in humans does not result in a fitness cost for invasive disease.

Mapping the Complex Genetic Landscape of Human Neurons.

When somatic cells acquire complex karyotypes, they are removed by the immune system. Mutant somatic cells that evade immune surveillance can lead to cancer. Neurons with complex karyotypes arise during neurotypical brain development, but neurons are almost never the origin of brain cancers. Instead, somatic mutations in neurons can bring about neurodevelopmental disorders, and contribute to the polygenic landscape of neuropsychiatric and neurodegenerative disease. A subset of human neurons harbors idiosyncratic copy number variants (CNVs, "CNV neurons"), but previous analyses of CNV neurons have been limited by relatively small sample sizes. Here, we developed an allele-based validation approach, SCOVAL, to corroborate or reject read-depth based CNV calls in single human neurons. We applied this approach to 2,125 frontal cortical neurons from a neurotypical human brain. This approach identified 226 CNV neurons, as well as a class of CNV neurons with complex karyotypes containing whole or substantial losses on multiple chromosomes. Moreover, we found that CNV location appears to be nonrandom. Recurrent regions of neuronal genome rearrangement contained fewer, but longer, genes.

Somatic nuclear mitochondrial DNA insertions are prevalent in the human brain and accumulate over time in fibroblasts.

The transfer of mitochondrial DNA into the nuclear genomes of eukaryotes (Numts) has been linked to lifespan in non-human species 1-3 and recently demonstrated to occur in rare instances from one human generation to the next 4. Here we investigated numtogenesis dynamics in humans in two ways. First, we quantified Numts in 1,187 post-mortem brain and blood samples from different individuals. Compared to circulating immune cells (n=389), post-mitotic brain tissue (n=798) contained more Numts, consistent with their potential somatic accumulation. Within brain samples we observed a 5.5-fold enrichment of somatic Numt insertions in the dorsolateral prefrontal cortex compared to cerebellum samples, suggesting that brain Numts arose spontaneously during development or across the lifespan. Moreover, more brain Numts was linked to earlier mortality. The brains of individuals with no cognitive impairment who died at younger ages carried approximately 2 more Numts per decade of life lost than those who lived longer. Second, we tested the dynamic transfer of Numts using a repeated-measures WGS design in a human fibroblast model that recapitulates several molecular hallmarks of aging 5. These longitudinal experiments revealed a gradual accumulation of one Numt every ~13 days. Numtogenesis was independent of large-scale genomic instability and unlikely driven cell clonality. Targeted pharmacological perturbations including chronic glucocorticoid signaling or impairing mitochondrial oxidative phosphorylation (OxPhos) only modestly increased the rate of numtogenesis, whereas patient-derived SURF1-mutant cells exhibiting mtDNA instability accumulated Numts 4.7-fold faster than healthy donors. Combined, our data document spontaneous numtogenesis in human cells and demonstrate an association between brain cortical somatic Numts and human lifespan. These findings open the possibility that mito-nuclear horizontal gene transfer among human post-mitotic tissues produce functionally-relevant human Numts over timescales shorter than previously assumed.

Analysis of Human Papilloma Virus Content and Integration in Mucoepidermoid Carcinoma.

Mucoepidermoid Carcinomas (MEC) represent the most common malignancies of salivary glands. Approximately 50% of all MEC cases are known to harbor CRTC1/3-MAML2 gene fusions, but the additional molecular drivers remain largely uncharacterized. Here, we sought to resolve controversy around the role of human papillomavirus (HPV) as a potential driver of mucoepidermoid carcinoma. Bioinformatics analysis was performed on 48 MEC transcriptomes. Subsequent targeted capture DNA sequencing was used to annotate HPV content and integration status in the host genome. HPV of any type was only identified in 1/48 (2%) of the MEC transcriptomes analyzed. Importantly, the one HPV16+ tumor expressed high levels of p16, had high expression of HPV16 oncogenes E6 and E7, and displayed a complex integration pattern that included breakpoints into 13 host genes including PIK3AP1, HIPI, OLFM4,SIRT1, ARAP2, TMEM161B-AS1, and EPS15L1 as well as 9 non-genic regions. In this cohort, HPV is a rare driver of MEC but may have a substantial etiologic role in cases that harbor the virus. Genetic mechanisms of host genome integration are similar to those observed in other head and neck cancers.

Early HPV ctDNA Kinetics and Imaging Biomarkers Predict Therapeutic Response in p16+ Oropharyngeal Squamous Cell Carcinoma.

PURPOSE: In locally advanced p16+ oropharyngeal squamous cell carcinoma (OPSCC), (i) to investigate kinetics of human papillomavirus (HPV) circulating tumor DNA (ctDNA) and association with tumor progression after chemoradiation, and (ii) to compare the predictive value of ctDNA to imaging biomarkers of MRI and FDG-PET. EXPERIMENTAL DESIGN: Serial blood samples were collected from patients with AJCC8 stage III OPSCC (n = 34) enrolled on a randomized trial: pretreatment; during chemoradiation at weeks 2, 4, and 7; and posttreatment. All patients also had dynamic-contrast-enhanced and diffusion-weighted MRI, as well as FDG-PET scans pre-chemoradiation and week 2 during chemoradiation. ctDNA values were analyzed for prediction of freedom from progression (FFP), and correlations with aggressive tumor subvolumes with low blood volume (TVLBV) and low apparent diffusion coefficient (TVLADC), and metabolic tumor volume (MTV) using Cox proportional hazards model and Spearman rank correlation. RESULTS: Low pretreatment ctDNA and an early increase in ctDNA at week 2 compared with baseline were significantly associated with superior FFP (P < 0.02 and P < 0.05, respectively). At week 4 or 7, neither ctDNA counts nor clearance were significantly predictive of progression (P = 0.8). Pretreatment ctDNA values were significantly correlated with nodal TVLBV, TVLADC, and MTV pre-chemoradiation (P < 0.03), while the ctDNA values at week 2 were correlated with these imaging metrics in primary tumor. Multivariate analysis showed that ctDNA and the imaging metrics performed comparably to predict FFP. CONCLUSIONS: Early ctDNA kinetics during definitive chemoradiation may predict therapy response in stage III OPSCC.

SquiggleNet: real-time, direct classification of nanopore signals.

We present SquiggleNet, the first deep-learning model that can classify nanopore reads directly from their electrical signals. SquiggleNet operates faster than DNA passes through the pore, allowing real-time classification and read ejection. Using 1 s of sequencing data, the classifier achieves significantly higher accuracy than base calling followed by sequence alignment. Our approach is also faster and requires an order of magnitude less memory than alignment-based approaches. SquiggleNet distinguished human from bacterial DNA with over 90% accuracy, generalized to unseen bacterial species in a human respiratory meta genome sample, and accurately classified sequences containing human long interspersed repeat elements.

SearcHPV: A novel approach to identify and assemble human papillomavirus-host genomic integration events in cancer.

BACKGROUND: Human papillomavirus (HPV) is a well-established driver of malignant transformation at a number of sites, including head and neck, cervical, vulvar, anorectal, and penile squamous cell carcinomas; however, the impact of HPV integration into the host human genome on this process remains largely unresolved. This is due to the technical challenge of identifying HPV integration sites, which includes limitations of existing informatics approaches to discovering viral-host breakpoints from low-read-coverage sequencing data. METHODS: To overcome this limitation, the authors developed SearcHPV, a new HPV detection pipeline based on targeted capture technology, and applied the algorithm to targeted capture data. They performed an integrated analysis of SearcHPV-defined breakpoints with genome-wide linked-read sequencing to identify potential HPV-related structural variations. RESULTS: Through an analysis of HPV+ models, the authors showed that SearcHPV detected HPV-host integration sites with a higher sensitivity and specificity than 2 other commonly used HPV detection callers. SearcHPV uncovered HPV integration sites adjacent to known cancer-related genes, including TP63, MYC, and TRAF2, and near regions of large structural variation. The authors further validated the junction contig assembly feature of SearcHPV, which helped to accurately identify viral-host junction breakpoint sequences. They found that viral integration occurred through a variety of DNA repair mechanisms, including nonhomologous end joining, alternative end joining, and microhomology-mediated repair. CONCLUSIONS: In summary, SearcHPV is a new optimized tool for the accurate detection of HPV-human integration sites from targeted capture DNA sequencing data.

Cas9 targeted enrichment of mobile elements using nanopore sequencing.

Mobile element insertions (MEIs) are repetitive genomic sequences that contribute to genetic variation and can lead to genetic disorders. Targeted and whole-genome approaches using short-read sequencing have been developed to identify reference and non-reference MEIs; however, the read length hampers detection of these elements in complex genomic regions. Here, we pair Cas9-targeted nanopore sequencing with computational methodologies to capture active MEIs in human genomes. We demonstrate parallel enrichment for distinct classes of MEIs, averaging 44% of reads on-targeted signals and exhibiting a 13.4-54x enrichment over whole-genome approaches. We show an individual flow cell can recover most MEIs (97% L1Hs, 93% AluYb, 51% AluYa, 99% SVA_F, and 65% SVA_E). We identify seventeen non-reference MEIs in GM12878 overlooked by modern, long-read analysis pipelines, primarily in repetitive genomic regions. This work introduces the utility of nanopore sequencing for MEI enrichment and lays the foundation for rapid discovery of elusive, repetitive genetic elements.

Expectations and blind spots for structural variation detection from long-read assemblies and short-read genome sequencing technologies.

Virtually all genome sequencing efforts in national biobanks, complex and Mendelian disease programs, and medical genetic initiatives are reliant upon short-read whole-genome sequencing (srWGS), which presents challenges for the detection of structural variants (SVs) relative to emerging long-read WGS (lrWGS) technologies. Given this ubiquity of srWGS in large-scale genomics initiatives, we sought to establish expectations for routine SV detection from this data type by comparison with lrWGS assembly, as well as to quantify the genomic properties and added value of SVs uniquely accessible to each technology. Analyses from the Human Genome Structural Variation Consortium (HGSVC) of three families captured ~11,000 SVs per genome from srWGS and ~25,000 SVs per genome from lrWGS assembly. Detection power and precision for SV discovery varied dramatically by genomic context and variant class: 9.7% of the current GRCh38 reference is defined by segmental duplication (SD) and simple repeat (SR), yet 91.4% of deletions that were specifically discovered by lrWGS localized to these regions. Across the remaining 90.3% of reference sequence, we observed extremely high (93.8%) concordance between technologies for deletions in these datasets. In contrast, lrWGS was superior for detection of insertions across all genomic contexts. Given that non-SD/SR sequences encompass 95.9% of currently annotated disease-associated exons, improved sensitivity from lrWGS to discover novel pathogenic deletions in these currently interpretable genomic regions is likely to be incremental. However, these analyses highlight the considerable added value of assembly-based lrWGS to create new catalogs of insertions and transposable elements, as well as disease-associated repeat expansions in genomic sequences that were previously recalcitrant to routine assessment.

Comprehensive identification of somatic nucleotide variants in human brain tissue.

BACKGROUND: Post-zygotic mutations incurred during DNA replication, DNA repair, and other cellular processes lead to somatic mosaicism. Somatic mosaicism is an established cause of various diseases, including cancers. However, detecting mosaic variants in DNA from non-cancerous somatic tissues poses significant challenges, particularly if the variants only are present in a small fraction of cells. RESULTS: Here, the Brain Somatic Mosaicism Network conducts a coordinated, multi-institutional study to examine the ability of existing methods to detect simulated somatic single-nucleotide variants (SNVs) in DNA mixing experiments, generate multiple replicates of whole-genome sequencing data from the dorsolateral prefrontal cortex, other brain regions, dura mater, and dural fibroblasts of a single neurotypical individual, devise strategies to discover somatic SNVs, and apply various approaches to validate somatic SNVs. These efforts lead to the identification of 43 bona fide somatic SNVs that range in variant allele fractions from ~ 0.005 to ~ 0.28. Guided by these results, we devise best practices for calling mosaic SNVs from 250× whole-genome sequencing data in the accessible portion of the human genome that achieve 90% specificity and sensitivity. Finally, we demonstrate that analysis of multiple bulk DNA samples from a single individual allows the reconstruction of early developmental cell lineage trees. CONCLUSIONS: This study provides a unified set of best practices to detect somatic SNVs in non-cancerous tissues. The data and methods are freely available to the scientific community and should serve as a guide to assess the contributions of somatic SNVs to neuropsychiatric diseases.

Genome diversity in Ukraine.

BACKGROUND: The main goal of this collaborative effort is to provide genome-wide data for the previously underrepresented population in Eastern Europe, and to provide cross-validation of the data from genome sequences and genotypes of the same individuals acquired by different technologies. We collected 97 genome-grade DNA samples from consented individuals representing major regions of Ukraine that were consented for public data release. BGISEQ-500 sequence data and genotypes by an Illumina GWAS chip were cross-validated on multiple samples and additionally referenced to 1 sample that has been resequenced by Illumina NovaSeq6000 S4 at high coverage. RESULTS: The genome data have been searched for genomic variation represented in this population, and a number of variants have been reported: large structural variants, indels, copy number variations, single-nucletide polymorphisms, and microsatellites. To our knowledge, this study provides the largest to-date survey of genetic variation in Ukraine, creating a public reference resource aiming to provide data for medical research in a large understudied population. CONCLUSIONS: Our results indicate that the genetic diversity of the Ukrainian population is uniquely shaped by evolutionary and demographic forces and cannot be ignored in future genetic and biomedical studies. These data will contribute a wealth of new information bringing forth a wealth of novel, endemic and medically related alleles.

A comparison of the psychological outcomes of self-reported and agency-notified child abuse in a population-based birth cohort at 30-year-follow-up.

BACKGROUND: Retrospective studies show a strong association between child abuse and subsequent psychiatric morbidity. Prospective studies are less common with conflicting data in young adults. We assessed the effect on psychological outcomes at 30-year follow-up of prospective agency notifications compared to retrospective self-reports of childhood maltreatment in the same birth cohort. METHODS: We used data on self-reported child abuse from 2425 young adults who completed the Child Trauma Questionnaire (CTQ) at the 30-year follow-up of a prospective birth cohort study commenced in 1981. These were linked to notifications of maltreatment to child protection agencies. The outcomes were DSM-IV diagnoses from the Composite International Diagnostic Interview-Auto version (CIDI). RESULTS: A quarter of the sample (n=600) self-reported maltreatment of any type, 326 (13.4%) rating this as severe. The most common maltreatment type was neglect (n=382), followed by emotional (n=225), sexual (n=198) and physical abuse (n=197). On adjusted analyses, there were differences between agency- and self-reported maltreatment in the psychological effects on 30-year-olds. In the case of depression, and especially post-traumatic disorder, there were significant associations irrespective of reporting sources. In the case of anxiety, there was a strong association with all forms of self-reported maltreatment. However, agency-reported cases were only significantly associated with emotional abuse in the previous 30-days. LIMITATIONS: The attrition rate from baseline to follow-up and the possibility of insufficient power to detect statistical significance in agency-reported cases CONCLUSIONS: The association between child maltreatment and psychiatric symptoms may vary by diagnosis and reporting source. Each source possibly captures different populations.

Association of CNVs with methylation variation.

Germline copy number variants (CNVs) and single-nucleotide polymorphisms (SNPs) form the basis of inter-individual genetic variation. Although the phenotypic effects of SNPs have been extensively investigated, the effects of CNVs is relatively less understood. To better characterize mechanisms by which CNVs affect cellular phenotype, we tested their association with variable CpG methylation in a genome-wide manner. Using paired CNV and methylation data from the 1000 genomes and HapMap projects, we identified genome-wide associations by methylation quantitative trait locus (mQTL) analysis. We found individual CNVs being associated with methylation of multiple CpGs and vice versa. CNV-associated methylation changes were correlated with gene expression. CNV-mQTLs were enriched for regulatory regions, transcription factor-binding sites (TFBSs), and were involved in long-range physical interactions with associated CpGs. Some CNV-mQTLs were associated with methylation of imprinted genes. Several CNV-mQTLs and/or associated genes were among those previously reported by genome-wide association studies (GWASs). We demonstrate that germline CNVs in the genome are associated with CpG methylation. Our findings suggest that structural variation together with methylation may affect cellular phenotype.

Long-term Cognitive, Psychological, and Health Outcomes Associated With Child Abuse and Neglect.

Potential long-lasting adverse effects of child maltreatment have been widely reported, although little is known about the distinctive long-term impact of differing types of maltreatment. Our objective for this special article is to integrate findings from the Mater-University of Queensland Study of Pregnancy, a longitudinal prenatal cohort study spanning 2 decades. We compare and contrast the associations of specific types of maltreatment with long-term cognitive, psychological, addiction, sexual health, and physical health outcomes assessed in up to 5200 offspring at 14 and/or 21 years of age. Overall, psychological maltreatment (emotional abuse and/or neglect) was associated with the greatest number of adverse outcomes in almost all areas of assessment. Sexual abuse was associated with early sexual debut and youth pregnancy, attention problems, posttraumatic stress disorder symptoms, and depression, although associations were not specific for sexual abuse. Physical abuse was associated with externalizing behavior problems, delinquency, and drug abuse. Neglect, but not emotional abuse, was associated with having multiple sexual partners, cannabis abuse and/or dependence, and experiencing visual hallucinations. Emotional abuse, but not neglect, revealed increased odds for psychosis, injecting-drug use, experiencing harassment later in life, pregnancy miscarriage, and reporting asthma symptoms. Significant cognitive delays and educational failure were seen for both abuse and neglect during adolescence and adulthood. In conclusion, child maltreatment, particularly emotional abuse and neglect, is associated with a wide range of long-term adverse health and developmental outcomes. A renewed focus on prevention and early intervention strategies, especially related to psychological maltreatment, will be required to address these challenges in the future.

Author Correction: A robust benchmark for detection of germline large deletions and insertions.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A robust benchmark for detection of germline large deletions and insertions.

New technologies and analysis methods are enabling genomic structural variants (SVs) to be detected with ever-increasing accuracy, resolution and comprehensiveness. To help translate these methods to routine research and clinical practice, we developed a sequence-resolved benchmark set for identification of both false-negative and false-positive germline large insertions and deletions. To create this benchmark for a broadly consented son in a Personal Genome Project trio with broadly available cells and DNA, the Genome in a Bottle Consortium integrated 19 sequence-resolved variant calling methods from diverse technologies. The final benchmark set contains 12,745 isolated, sequence-resolved insertion (7,281) and deletion (5,464) calls ≥50 base pairs (bp). The Tier 1 benchmark regions, for which any extra calls are putative false positives, cover 2.51 Gbp and 5,262 insertions and 4,095 deletions supported by ≥1 diploid assembly. We demonstrate that the benchmark set reliably identifies false negatives and false positives in high-quality SV callsets from short-, linked- and long-read sequencing and optical mapping.

Biological diversity in an Islamic archaeological population: A radiogenic strontium isotope and craniometric analysis of affinity in Ottoman Romania.

OBJECTIVES: Written accounts, as well as a previous craniometric study, indicate that migrations of non-Europeans and conversions of Europeans to Islam define Ottoman communities in Early Modern Europe. What is less clear are the roles of migration and admixture in generating intra-communal variation. This study combines craniometric with strontium isotope data to compare the cranial affinities of locally born and immigrant individuals. We predict that locally born individuals are more likely than non-locals to show evidence of admixture. MATERIALS AND METHODS: Radiogenic strontium isotope data for 21 Ottomans were compared against archaeological faunal values. Sixteen individuals with intact crania were also measured and compared against two comparative source populations from Anatolia and Europe. Discriminant function analysis assigned unclassified Ottoans to either comparative group based on typicality probabilities, with potential admixture established via intermediate morphology between the two source populations. RESULTS: Strontium isotope values revealed relatively high proportions of non-locals, consistent with high mobility documented historically. The sexes differed, with more males classifying as "typically Anatolian" than females. Locals and non-locals also had different cranial affinity patterns, with most classifying either as "typically Anatolian" or "typically European." Contrary to expectation, none of the locals were identified as intermediate, suggesting admixture rates were relatively low. CONCLUSIONS: Consistent with historical records, the results revealed high levels of extra-regional migration, with most individuals identifiable as either typically Anatolian or European. Moreover, locals and non-locals differed craniometrically, with no signs of admixture between Anatolian migrants and European converts in locals. This suggests intra-communal divisions were maintained.

Characterization of nuclear mitochondrial insertions in the whole genomes of primates.

The transfer and integration of whole and partial mitochondrial genomes into the nuclear genomes of eukaryotes is an ongoing process that has facilitated the transfer of genes and contributed to the evolution of various cellular pathways. Many previous studies have explored the impact of these insertions, referred to as NumtS, but have focused primarily on older events that have become fixed and are therefore present in all individual genomes for a given species. We previously developed an approach to identify novel Numt polymorphisms from next-generation sequence data and applied it to thousands of human genomes. Here, we extend this analysis to 79 individuals of other great ape species including chimpanzee, bonobo, gorilla, orang-utan and also an old world monkey, macaque. We show that recent Numt insertions are prevalent in each species though at different apparent rates, with chimpanzees exhibiting a significant increase in both polymorphic and fixed Numt sequences as compared to other great apes. We further assessed positional effects in each species in terms of evolutionary time and rate of insertion and identified putative hotspots on chromosome 5 for Numt integration, providing insight into both recent polymorphic and older fixed reference NumtS in great apes in comparison to human events.